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Experimental  demonstration  of  modern  modal  control  of  a 
flexible  clamped-clamped  beam  in  bending  was  accomplished. 
Pole  placement  techniques  in  conjunction  with  a  Leuenberger 
Observer  were  used  for  controller  development.  Available 
microcomputer  computational  capabilities  drove  the  reduced 
order  model  to  a  two  mode  representation.  The  controllers, 
designed  with  varying  degrees  of  robustness,  demonstrated  the 
limitations  of  reduced  order  controllers  in  even  the  simplest 
of  systems.  Sensor  and  actuator  position  are  demonstrated  to 
be  of  paramount  importance  in  the  modern  controller  LQ  design. 


iTRQDUCTION 


The  space  shuttle  is  operational.  The  USSR  has  a 
permanent  presence  in  space.  Private  industry  is  well  on  its 
way  to  having  its  own  satellite  launching  facility.  These 
recent  developments  can  only  mean  one  thing;  trends  towards 
larger,  and  more  complex  space  structures  are  quickly  becoming 
a  reality.  Demands  in  pointing,  jitter,  and  surface 
tolerances  for  these  satellites  will  require  more  robust  and 
adaptive  controllers  due  to  inaccuracies  in  the  mathematical 
models  and  changing  mass  distribution.  Designing  a  controller 
tor  these  satellites  amounts  to  attempting  to  control  the 
minimum  number  of  structural  modes  in  order  to  maintain 
desired  shape  and  attitude  control  requirements. 
Unfortunately,  this  reduced  order  controller  is  affected  by 
the  structure's  real  modes  not  included  in  its  model.  Their 
effect  on  control  inputs  is  called  control  spillover  and  their 
effect  on  sensor  outputs  is  known  as  observation  spillover. 
Both  can  cause  structural  instabilities  which  exceed  the 
design  requirements. 

Calico  and  Janiszewski  (Ref  1)  demonstrated  that 
elimination  of  either  observation  or  control  spillover  will 
ensure  system  stability.  The  number  of  modes  whose 
observation  spillover  can  be  eliminated  is  one  less  than  the 
number  of  sensors,  while  the  number  of  modes  whose  control 
spillover  can  be  eliminated  is  one  less  than  the  number  of 
actuators.  It  becomes  obvious  at  this  point  that  the  most 


oractical  implementation  will  be  elimination  of  observation 


spillover  because  of  the  physical  limitations  of  implementing 
a  myriad  of  actuators. 

Computer  simulation  of  this  phenomena  has  been 
demonstrated  repeatedly  (Ref  1),  and  controller  design  for 
large  space  structures  has  progressed  rapidly  in  recent  years. 
These  developments  are  now  available  for  demonstration  of 
control  on  physical  structures.  Schaechter  (Ref  2)  has 
demonstrated  active  shape  and  dynamic  control  by  using  various 
control  law  design  approaches,  including  adaptive  control,  and 
demonstrated  various  microprocessor  control  capabilities  on  a 
very  flexible  clamped-free  12.5  foot  beam.  However,  the 
requirements  for  sensor  and  actuator  sensitivity,  as  well  as 
computational  capability  preclude  the  use  of  such  a 
sophisticated  set-up  to  demonstrate  active  modal  control. 

The  approach  used  in  this  experiment  closely  parallels  a 
control  demonstration  by  Herrick  (Ref  3)  on  a  flexible,  5  foot 
clamped-clamped  beam.  This  simple  case  is  a  reasonable 
starting  point  for  demonstration  of  active  modal  control. 
Even  for  this  simple  case,  computer  hardware  available  impacts 
the  structure  and  complexity  of  the  test  significantly.  The 
result  is  a  control  design  based  on  a  significantly  reduced 
order  model  and  an  observer  based  on  only  a  single 
measurement. 

The  objectives  of  Liiis  experiment  are  to;  1) 
demonstrate  active  modal  control;  2)  demonstrate  effects  of 
observation  spillover  of  residual  modes  not  modeled  in  the 
controller;  and  3)  provide  a  basis  tor  further  modal  control 
testing  at  the  Air  Force  Institute  of  Technology  (AFIT). 


The  laboratory  developed  to  demonstrate  active  modal 
control  consist  of  the  test  bed  and  beam,  modal  identification 
equipment,  and  control  implementation  equipment. 

The  test  bed  is  very  massive  (over  one  ton)  with  rigid 
beam  clamps  and  rigid  actuator  mounts.  The  beam  is  mounted  on 
edge,  and  the  actuators  provide  force  perpendicular  to  it  in 
the  horizontal  direction.  It  was  desireable  to  keep  the  first 
mode  of  vibration  above  10  hz  since  the  lightest  accelerometer 
available  (  2  gm  )  was  only  calibrated  down  to  that  frequency. 
The  beam  chosen  therefore  was  a  61.5“  x  1"  x  .25“  beam  of 
2024-T-4  aluminum,  clamped  at  both  ends.  The  accelerometer 
is  the  piezoelectric  type  with  a  preamplifier  within,  and  a 
selectable  gain  amplifier  also  provided.  The  actuator  is  a 
typical  classroom  shaker  employing  a  voice  coil  system  and  has 
a  sizeable  shaker  core  mass(341gm).  The  test  bed  is  versatile 
in  its  ability  to  reconfigure  present  equipment  and  can  also 
accommodate  mounting  of  larger  and  more  complex  test 
structures.  The  details  of  the  equipment  set-up  are  located 
in  Appendix  1. 

The  equipment  used  for  modal  parameter  identification 
consists  of  the  apparatus  used  for  sine  dwell  identification 
and  the  apparatus  used  for  spectrum  analysis.  Sine  dwell 
equipment  includes  a  digital  frequency  counter,  two  root  mean 
square  voltmeters,  a  storage  oscilloscope,  and  a  sine  function 
generator.  The  spectrum  analysis  utilizes  a  band  limited 
white  noise  generat  r,  a  signal  correlator,  and  a  spectrum 


display.  The  spectrum  display  includes  storage  display  of 
power  spectral  density  and  transfer  to  plotter  for  hardcopy 
results. 

The  sensor  and  actuator  just  described  were  utilized  in 


the  system  to  implement  the  control  algorithm  shown  in  Fig  1. 


Figure  1.  Laboratory  Control  Components. 

The  accelerometer  is  calibrated  periodically,  including 
its  wires, and  the  tiny  operational  amplifier  inside  also 
helps  guarantee  accurate  sensitivity  and  reduced  noise.  A 
vibration  meter  specifically  designed  for  vibration  testing 
incorporates  a  5  Hz  high  pass  filter  and  outputs  the  desired 
displacement  signal. 

The  A/D  and  D/A  converters,  memory,  and  single  board 
computer  were  housed  together  and  access  was  obtained  by  a 
video  terminal.  Floppy  disk  program  storage  was  essential  to 


both  software  development,  implementation,  and  data  analysis. 
A  simple  set  of  potentiometers  attenuated  the  desired 


force  signal  so  that  a  summer  could  compare  it  to  the  force 
output  signal  of  the  impedance  head.  The  pot  settings  could 
be  very  accurately  controlled  using  the  analog  computer's 
digital  readout  in  the  pot  set  mode.  This  became  helpful  when 
compensating  for  the  input  scale  factor  in  the  digital  control 
implementation,  in  addition,  the  accelerometer  within  the 
impedance  head  could  have  provided  negative  damping  to  more 
closely  approximate  the  true  beam's  dynamics.  However, 
measured  beam  damping  closely  matched  the  light  predicted 
damping  in  modal  identification  testing  and  was  therefore  not 


used. 
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EXPERIMENTAL  PROCEDURE 


The  experimentation  consists  of  modal  parameter 
identification,  computer  software  development,  and  finally  the 
control  experiment. 

The  beam’s  modal  parameters  of  concern  are  natural 
frequencies  and  damping.  Two  methods  were  used  to  identity 
these  parameters;  swept  sine  and  spectrum  analysis. 

Swept  Sine,  using  the  set-up  decribed  in  detail  in 
Appendix  1,  a  single  actuator  was  driven  using  a  constant  root 
mean  square  voltage.  Then  the  magnitude  of  the  transfer 
function  between  the  force  in  and  the  acceleration  out  was 
calculated  for  a  range  of  frequencies  near  each  natural 
frequency.  The  frequency  was  determined  accurately  from  the 
digitally  displayed  signal  period.  Then  the  magnitude  of  the 
transfer  function  was  plotted  versus  frequency,  and  the  peaks 
were  used  to  identify  the  first  tour  natural  frequencies  (to-) 
of  vibration.  The  half  power  method  was  then  applied  to 
determine  the  damping  ratio  ( ). 

Spectral  Analysis.  The  shaker  was  driven  by  a  bandlimited 
white  gaussian  noise  generator  which  inputs  equal  power  to  all 
frequencies  within  its  bandwidth.  Then  a  number  of  samples  (on 
the  order  of  2expl5)  of  the  accelerometer  output  were 
correlated  with  a  controllable  period,  T,  between  samples.  The 
spectrum  analyzer  takes  the  Fourier  Transform  of  this 
autocorrelation  and  displays  the  accelerometer  signal's  power 
spectral  density.  The  sensitivity  of  the  frequency  scale,  Af, 
displayed  on  the  horizontal  axis  is  a  function  of  the  sampling 
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period,  T: 


At  =  1/(20T)  (Hz/cm)  (1) 

Typical  values  for,  T,  are  lOms,  3.33msec,  and  1msec  which 
allow  frequency  ranges  of  0-50Hz,  O-lSOHz,  and  0-500Hz.  Care 
must  be  taken  that  the  bandwidth  of  noise  input  does  not 
excite  beam  natural  frequencies  higher  than  the  sampling 
ability  of  the  correlator.  This  would  manifest  itself  in  the 
form  of  false  spectral  content  in  the  frequencies  within  the 
bandwidth  of  the  spectrum  analyzer(aliasing). 

The  natural  frequencies  and  damping  were  obtained  from  the 
plots  of  power  spectral  density  in  a  similar  manner  as  with 
swept  sine  except  that  now  the  PSD  of  the  output  is  plotted 
versus  frequency  where; 

2 

PSD  y(t)  =  |H(jiO)|  X  PSD  f(t)  (2) 

In  Eq  2,  f(t)  is  the  force  input  signal,  y(t)  is  the 
accelerometer  output  signal,  and  (H(ju>)|  is  the  magnitude  of 
the  transfer  function. 

Background  noise  was  also  examined  in  this  manner  to 
determine  the  significance  of  its  effect  on  the  system.  A 
coherence  and  Modal  Assurance  Criteria  analysis  was  not 
required  because  the  beam  natural  frequencies  were  expected  to 
be  widely  separated  and  lightly  damped. 

Software  Development.  Computer  software  was  developed  to 
improve  the  speed  of  the  single  board  computer  (on  the  S-100 
bus)  from  2.5  MHz  to  5.0  MHz  operation.  Logic  was  then 
implemented  in  machine  code  to  obtain  accurate  time 
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measurement  trom  the  computer's  5  MHz  system  clock.  Next,  the 
machine  language  program  for  operation  of  the  analog  to 
digital  board  and  the  digital  to  analog  board  was  developed. 
The  boards  were  calibrated  and  then  checked  using  simple 
multipl  channel  input-output  routines.  A  first  attempt  at  a 
control  algorithm  was  written  in  Fortran,  but  was  found  to  be 
much  to  slow  to  obtain  the  desired  sampling  rate  of  200  Hz. 
This  rate  was  chosen  to  obtain  approximately  six  samples  per 
cycle  of  the  second  vibrational  mode  of  the  beam. 

A  number  of  techniques  should  be  noted  here  which  can 
quickly  reduce  the  computation  time  for  real  time  data 
processing  beyond  hardware  improvements: 

1)  Use  integer  arithmetic  where  it  is  possible  to  do 
so  and  still  obtain  the  desired  accuracy, 

2)  Avoid  inefficient  computation  logic  and  divides, 

3)  Use  the  assembly  language  version  of  the  higher 
order  language  program  to  identify  ineffective 
and  redundant  commands  created  by  the  compiler. 

4)  Avoid  the  excessive  use  of  subroutine  calls,  data 
storage  look-ups,  and  unnecessary  16-bit  loads. 

After  optimization,  however,  the  control  algorithm 
became  specific  to  the  number  of  modes  controlled,  the  number 
of  sensors,  the  number  of  actuators,  and  the  scaling  factors 
required  to  avoid  signed  16-bit  overflow.  Future 
experimentation  will  require  scheduled  events  to  occur  in  an 
intricate  manner.  A  series  of  programs  for  control 
implementation  were  developed  for  this  purpose  and  can  be 
found,  along  with  a  number  of  test  programs,  in  the  laboratory 


a 


as  documented  in  Appendix  8,  The  control  implementation 
program  and  its  subroutines  are  contained  in  Appendix  5. 

The  Control  Experiment.  The  control  was  applied  and 
evaluated  for  two  cases:  initial  condition  input  and 
continuous  vibration  input.  The  initial  condition  applied  was 
a  step  input  using  a  static  deflection  of  the  beam  center  of 
approximately  one  third  of  an  inch.  A  digitized  temporal 
estimate  of  the  first  two  modal  amplitudes  was  the  analysis 
tool  used  to  compare  modal  response  parameters  without 
control,  with  only  the  force  loop  around  the  shaker  active, 
and  with  control  applied  (for  details,  see  Appendix  1). 

The  continuous  vibration  applied  was  white  gaussian  noise 
limited  to  ,  SOHz,  150Hz,  and  500Hz  at  RMS  voltage  levels  of 
,1,  .2,  and  .3  volts.  This  was  applied  through  a  second 
shaker  positioned  at  the  end  opposite  to  the  control  shaker 
and  equidistant  from  the  beam  center.  Incorporating  the 
results  of  previous  analysis  by  Hungerford  (Ref  4),  care  was 
taken  not  to  place  either  actuator  near  a  node  of  the  first 
four  mode  shapes  to  avoid  diminished  control  authority.  Also 
to  avoid  large  actuator  deflections  which  might  cause 
nonlinearities,  the  actuators  were  not  centrally  located.  The 
sensor  was  also  located  away  from  the  nodes  to  avoid  a 
theoretically  infinite  response  time.  Controlled  modes  and 
observation  spillover  terms  from  the  third  and  fourth  modes 
were  then  examined  by  spectral  analysis  for  the  various  cases, 
and  compared  to  theoretical  calculations '  expected  root 
migration  due  to  the  control  applied. 
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THEORETICAL  DEVELOPMENT 

The  development  of  the  control  algorithnl,  its  elements, 
and  its  effects  on  other  system  components  must  begin  with  the 
test  beam.  The  controller  is  then  designed,  and  problems 
involving  its  implementation  are  also  addressed  here.  Finally 
the  equipment  sensitivities  are  calibrated  and  matched. 


The  Test  Beam.  Physical  properties  are  identified  below; 
Material:  2024  T-4  Aluminum 

Length;  61.5"  =  1  Nondimensi ional  Length  =  1  NL 

Actuator  Position;  x  =  4.9170/61.5  =  0.079956  NL 

a 

Sensor  Position;  x  =  48.468/61.5  =  0,788106  NL 

s 


Cross-Section:  1.0"  x  0.25" 

3 


Density;  yo  =  0,097  Ibm/in 
Modulus  of  Elasticity;  E  = 


6 

10.1  X  10 
3 


Area  Moment  of  Inertia;  I  *  bh  /12  = 
Mass  per  unit  length;  m  =  yoA/(12g) 


2 

Ibf/in 

4 

001302  in 


=  0.097(.25)/32.2(12) 

-5  2  2 

=  6.27588  X  10  Ibf-sec/in 


Natural  Modes.  It  was  first  necessary  to  determine  the 
natural  modes  and  frequencies  for  a  uniform  clamped-clamped 
beam  in  bending.  The  general  equation  of  motion  for  a  beam  in 
bending  is; 


»  X* 


[  El 


T7 


]  =  -  m 


^  t* 


(3) 
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where  the  form  of  v(x,t)  using  the  assumed  modes  technique  is: 


v(x,t)=  V  V  (x)v  (t) 
Z_i=l  i  i 


(4) 


For  a  uniform  beam  the  spatial  coordinates  v  (t)  satisfy: 

i 


V  ""  -  (^A/EI)V  =  0 
i  i  '  i 


(5) 


The  solution  to  this  homogeneous  equation  is: 


V  (x)=  C  sin(|3  x)  +  C  cos(A  x)  +  C  sinh(|3  x)  +  C  cosh(j3  x)  (6) 
ili2i3i4i 


2  4 

where  =  A/3  /El 

i  i 


(7) 


Boundary  conditions  for  the  clamped-c lamped  case  arei 


V  (0)  =  V  (0)  =  V  (L)  =  V  (L)  =  0 
i  i  i  i 


(8) 


Applying  the  boundary  conditions,  we  have: 


V  (0)  =  0  =  C  +  C 
i  2  4 


v'  (0)=  0  =  C  +  C 
i  3  4 

(9) 

V  (L)  =  0  =  C  sin(|3  L)+C  cos(Q  L)+C  sinh(fl  L)+C  cosh(A  L) 
i  112  ”i  3  ^  i  4  ”i 


9 

V  (L)=  0  =  [C  cos(i3  L)-C  sin()3  L)+C  cosh(A  L)+C  sinh(A  L)1 
i  11^12  3  *  i  4  *  i 


Combining  these  relationships  gives: 


[sin(|3  L)-sinh(^  L)]  [cos(^  L)-cosh(^  L)] 
i  i  i  i 

[cos(^  L)-cosh(|3  L)]  [-sin(|J  L)--sinh(/3  L)] 
i  i  i  i 


Zeroing  the  determinant  is  equivalent  to: 


cos(^  L)cosh(j3  L)  =  1  (11) 

'  i  i 


From  Eq  (7)  we  find  the  solutions  to  Eq  (11)  to  be: 

2  4*5 

lO  =  (P  L)  [EI/mL  ]  (12) 

i  i 


where 


(Q  L)  =  4.730 

1 

(/3  L)  = 
2 

7.853 

(a  L)  =  10.996 

3 

(A  L)  = 

4 

14.137 

and  for  larger  values  of  i: 

{p>  L) 

=  (2i+l)'Tt/2 

Assuming  lengths  will 

now  be  thought  of 

in  terms  of  beam 

lengths,  the  theoretical 

values  for  the 

beam's  natural 

frequences  are  for  L  =  1  NL: 


/3 

1 

=4.73: 

U)  =85.6278 
1 

rad/sec 

or  13.6281 

Hz 

(9 

2 

=7.853; 

U3  =236.028 
2 

rad/sec 

37.565 

Hz 

(3 

3 

=10.996: 

U)  =462.767 
3 

rad/sec 

73.652 

Hz 

/3 

=14.1372 

to  =764.922 

rad/sec 

121.740 

Hz 

4  4 
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Common  practice  is  to 


Normalizing  The  Eigenvectors, 
normalize  with  respect  to  the  mass.  This  requires: 


PL  2 

\  m  V  (x)  =  1  i=l,2,3, . . .n  (13) 

JO  i 


So, 


V  (x)=C  sin(A  x)-C  sinh(/3  x)+  Z  C  cos(A  x)-  Z  C  cosh(A  x)  (14) 
i  1  il  i  il  i  il  i 


where 


C  =  Z  C  (15) 

2  i  1 


Solving  for  C  and  C  from  Eq  (8)  and  Eq  (13)  (in  Appendix  2); 
1  2 

C* 


For 

i*l : 

C  =  .80452 

1 

and 

C  = 
2 

-.81885 

For 

i=2: 

C  =  .81926 

1 

and 

II 

CM 

u 

-.81862 

So  that 


V  (x)  =  .80452 [sin( 4 .730x)  -  sinh( 4 .730x) ] 

1 

+.81885[cosh(4.730x)  -  cos(4.730x)]  (17) 

V  (x)  »  .819261lsin(7.853x)  -  sinh(7 .853x) ] 

2 

+.81862[cosh(7.853x)  -  cos(7.853x)]  (18) 


LOG  Regulator  Development.  Under  the  assumptions  of 


linear  models,  quadratic  cost  criterion,  and  white  Gaussian 
noise  inputs,  the  design  of  the  optimal  controller  can  be 
completely  separated  from  the  design  of  the  estimator.  Using 


Hungerford's  development  of  the  beam  model  developed  by  Balas 
(Ret  6),  the  optimal  control  problem  can  be  expressed  by  the 
state  equations: 


Xc  = 

[Acl 

Xc  + 

[Be] 

U 

• 

Xr  = 

[Ar] 

Xr  + 

[Br] 

U 

where 


T 

• 

• 

XC  = 

(  V  V  V 

V  ) 

12  1 

2 

T 

• 

• 

xr  = 

(  V  V  V 

V  ) 

3  4  3 

4 

and  where 


(19) 


(20) 


[AC] 


I 

0  0  j  1  0 

0  0  I  0  1 

— I-— 

2  j 

-t)  j-2tu> 

•i  •  i 


i  =  1,2  (21) 


[Be]  =  (00  V  (X  )  V  (x  )  ] 

la  2  a 


(22) 


and  where 


[Ar] 


0 

0 

1  0 

0 

0 

_ J 

0  1 

2 

-U) 

-2tu) 

i 

i 

3,4  (23) 
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[Br]  =  [00  V  (x  )  V  (x  )  ] 

la  2  a 


(24) 


In  addition,  the  output  equations  are  given  by 


Y  -  [Cc]  Xc  +  [CrJ  Xr 


where 


[Cc]  =  (V  (x  )  V  (x  )  0  0] 
Is  2  s 


(26) 


[Cr]  =  [  V  (x  )  V  (x  )  0  0  ] 

3  s  4  s 


The  quadratic  cost  function  to  be  minimized  is; 


J  » 


fS 


oo  _  T  _  _  T  _ 

(X  lOJ  X  +  U  [R]  U  )  dt 

0 


The  local  matrix  Riccati  equation  solving  package  OPTCON 
(Ref  5)  solves  for  the  feedback  matrix  [G]  so  that  Eq  (28)  is 
minimized,  within  the  constraints  of  Eq  (27), by  the  control 
applied  U  : 


U  *  -  [G]  X 


where 


-1  T 

[G]  =  IR]  [B  ]  [Si 


[S]  is  the  solution  to  the  steady-state  Ricatti  equation: 


T  -IT 

[S]  [Ac]  +  [Ac]  [S]  -  [S]  [Be]  [R]  [Be]  [S]  +  [Q]  »  [0]  (31) 
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The  OPTCON  package  solves  the  Ricatti  equation  with 
the  appropriate  values  for  [Acl ,  [Be],  IQI ,  and  [R].  Four 
cases  were  considered  for  lOi  :  [Q]  =  20,  200,  500,  and  6000. 
For  now,  we  will  consider: 

[01=200[I]  and  [R]  =  [I]  (32) 

If  full  state  feedback  is  to  be  implemented,  a  state 
estimator  will  be  required  to  produce  an  estimate  of  Xc  from 
the  output  Y  .  The  deterministic  state  estimator  is  of  the 
form: 


A 

x(t)c=  [Ac]  x(t)c  +  [Bc]u(t)  +  [K]  [  Y(t) 


?(t) 


(33) 


X(t)c  is  the  state  vector  estimate,  y(t)  is  the  sensor  output, 
and  Y(t)  the  output  corresponding  to  the  estimated  state. 
U(t)  is  the  actuator  input.  The  optimum  observer  gain  matrix 
[K]  ,  can  be  solved  for  by  the  same  OPTCON  package  with  the 
following  substitutions: 


Replace 


[Ac] 

with 

[AC] 

[BC] 

with 

[Cc] 

The  output  is  the  matrix  of  estimator  gains,  [K] ,  transposed. 

The  resulting  regulator  in  Figure  2  is  reduced  to  a 
precomputable  set  of  gains  to  be  applied  to  the  output  Y(t)  to 
obtain  the  desired  input  to  the  shaker  U(t). 


Figure  2.  Control  Block  Diagram 
Therefore,  the  desired  control  to  apply  ist 

-1 

0(t)  -  IGJ {8{IJ-[AcJ~{BcJ {GI+IKJ (CcJl  (Kj  Y(t)  (34) 

The  stability  of  the  controller  is  therefore  governed  by 


the  eigenvalues  of  the  matrix  [A] .  Where  [A]  is  given  by 

[A]  =  [[AC]  +  tBc] [G]  -  lK][Cc]]  (35) 


The  computed  natural  frequencies , u>  ,  which  are  slightly 

i 

higher  than  the  measured,  are  used  for  the  formulation  of 


[Ac]  : 


[AC]  = 


0 

0 

1 

0 

0 

0 

0 

1 

-7332.12 

0 

-.8563 

0 

0 

-55709.0 

0 

-2.3603 

(36) 


[Be]  =[00  .10262  .25397  ] 

[Cc]  =  [  .307124  -.694532  0  0  ] 

[G]  =  [  -2.3698  7.2801  8.0764  7.6288  ] 

T 

[K]  =  [  11.6250  -11.145  -23.573  -19.485  ] 
and,  as  computed  (details  are  also  in  Appendix  3)  [A]  is: 


(37) 

(38) 

(39) 

(40) 


[A]  = 


-3.572 

8.077 

!  1.000 

0 

3.424 

-7.744 

I  0 

1 

1.000 

-7325. 

-15.62 

I-.0275 

.7829 

15.384 

-55720 

1  2.051 

1 

-.4228 

(41) 


[A]  has  eigenvalues:  -1.800  +/“  J  85.63 

-4.083  +/-  J  235.9 

The  system  eigenvalues  are  represented  by  the  eigenvalues 
of  [[Ac]  +[Bc][G]]  and,  as  computed  in  Appendix  3,  are: 


-.8426  +/-  85.624 

-2.149  +/-  236.02 
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h  -  - 

L' 


i 


Q 


As  developed  by  Calico  and  Janiszewski  in  Ref  1,  the 
overall  system  stability  with  residual  modes  present  is 
represented  in  Eq  (42),  Note  that  the  off-diagonal  terms 
represented  by  [K] [Cr]  from  Eq  (27),  are  due  to  observation 
spillover,  and  those  represented  by  [Br] [G] ,  from  Eq  (24), 
are  due  to  control  spillover.  The  actual  matrix  elements  and 
matrix  eigenvalues  for  various  IQ]  values  are  in  Appendix  3. 


Z 


[  [Ac]  +  [Bc]  [G] ] 
[0] 

[[Br] [G]] 


[[Be] [G] 1 
[(Ac]-tK] ICcl] 
t[Br] [GJ] 


toi 

t[K) [Cr]] 
[Ar] 


Z 


(42) 


and  the  augmented  state  vector,  Z  is: 


Z  = 


where 


e  (t)  =  V  (t) 
c 


V  (t) 

c 


(43) 


(44) 


For  the  sensor  and  actuator  locations  chosen,  [Cr]  and  [Br] 
are  given  by; 


[cr]  =  [  1.0172  -1.2294  0  0  ]  (45) 

[Br]  =[00  .4485  .6627  ]  (46) 


1 
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and  [Ar]  is 


[Ar]  = 


0 

0 

! 

0 

0 

0 

0 

1 

_i _ 

0 

0 

-214153 

0 

1 

1  -4. 

6277 

0 

0 

-585106 

1 

1 

0 

-7.6492 

The  equations  tor  the  estimator  and  controller  are  thus: 


^  A 

X  =  [A]  X  +  [K]Y 


U  =  [G]  X 


(48) 

(49) 


where  Y  is  the  measurement,  and  U  is  the  control  to  apply. 

Control  Implementation,  To  obtain  a  discete  ot 

Eq  (34),  we  tirst  integrate  both  sides: 


A  lA]  t  A  rt  lA]  (t-Tr) 

X(t)  =  e  X(0)  +  \  e  [K]Y(T:)  dX 


or,  as  derived  in  Appendix  4, 


A  [A]T  ^  -1  [A]T 

X(T)  =  e  X(0)  +  [A]  [e  -  [I]][KlY(0) 


The  result  of  expanding  the  exponential  into  a  series,  is: 


A  lAjT  A  -1  ([AJT) 

X(T)  =  e  X(0)  +  lA]  [[A]T  +  - 

21 


...)  [Kl  Y(0)  (52) 


In  general,  for  the  kth  sample  period, 


A  ([A]T)  A 

X(T  )  =  [[I]  +  [A]T  +  -  +  ...]  X(T  ) 

k  21  k-1 


(A]T  ([A)T) 

+  T[[I]  +  —  +  -  '  ■—  +  ...][K]y(T  ) 

21  31  k-1 


(53) 


The  discrete  control  to  be  implemented  by  the  computer  is: 


A  A 

X(T  )  =  [H(T)]  X(T  )  +  [F{T)]y(T  )  (54) 

k  k-1  k-1 


Or  simply, 

A'k  A  A 

X  =  [H]  X  +  (F]y  (55) 

i+1  i  i 


A 

U  =  IG]  X 
i+1  i+1 


where 

2  3 

(IA)T)  (IA]T) 

(H(T)]  =  [I]  +  IA)T  +  -  +  -  +  ... 

21  31 


(56) 


(57) 


and 


2 

lAjT  (lA)T) 

[F(T))  =  T[II]  +  -  +  -  +  ...]lK]  (58) 

21  31 


Using  a  sample  period  of  16  clock  ticks  (clicks)  or  about 

T  =  .008  sec 
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Due  to  the  limited  floating  point  capacity  of  the  single 
board  computer,  these  matrix  multiplications  must  be  carried 
out  in  integer  arithmetic  to  adequately  follow  the  arriving 
frequencies  at  the  sensor.  Note  that  [H]  contains  a  broad 
range  of  values  among  its  elements.  The  range  is  so  broad 
that  to  use  the  matricies  as  is  would  cause  the  computed 
values  to  exceed  the  maximum  value  allowed  in  signed  16-bit 
integer  arithmetic.  Thus  Eq  (55)  is  carefully  scaled: 


Where  ^  can  be  thought  of  as  a  new  "pseudostate"  vector. 
Now  Eqs  (55)  and  (56)  become: 
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•is^  .7513  .0243  .0073  .0002 

‘^^2  =  -.3220  .0001  .0040  .1^ 

10^^  -b3.38  -2.316  .7765  -.0014  10  ^ 

-4.493  -220.9  -.0051  -.2936  lo  ^ 


.in 

+  ^  2  Y  (62) 

lOf] 


f-2.3698  7.2801  8.0764  7.6288] 


2  = 


.7513  .0243  .7287 
.0108  -.3220  .0121 
5338  -.0232  .7765 
0449  -2.209  -.0051 


.0193 


.3964 


-.0014 


-.2936 


2  + 


.8205 


.431( 


Y  (64) 
.2688  i 

1.414 


f-.2370  .7280.  80.76  76.24] 
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To  make  these  equations  integer  valued,  multiply  [H]  and 

[F]  by  a  factor  of  128,  and  divide  ^  by  128  when 

i  +  1 

complete.  Thus  we  define  a  new  pseudostate  vector?,  ^  : 


96  3  93 


1  -41  2  51 


-68  -3  99  0 


-6  -283  -1  -38 


Now  the  original  pseudostate  vector  is: 


and  the  control  to  apply  over  the  next  interval  is: 


U  =  [0  1  81  76] 

i-H 
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Thus  the  control  calculations  reauce  to  a  fairly  simple 
set  of  multiplies  ana  aiviaes  in  the  actual  algorithm 
(Appendix  5). 

Sensitivity  Analysis.  Sample  calculations  indicate  that 
these  equations  will  not  exceed  the  maximum  integer  value  of 
2expl5  as  long  as  the  measurement  rangeis  constrained  to 
integers  with  Yi;  -100<Yi<+100 

The  analog  to  digital  converter  is  set  up  such  that  a 
linear  relationship  exists  over  the  range  Vin:  -10v<Vin<+10v 
where : 


Vin  =  Voltage  to  the  A/D  board 

Din  =  Digitized  value  of  input  voltage 


so  that 


Din  =  (Vin+10)  x  204.8 


Thus  tor: 


Vin  = 

-10 

=> 

Din  = 

0 

Vin  = 

0 

=> 

Din  = 

2048 

Vin  = 

+  10 

=> 

Din  = 

4096 

Now  defining  In  =  Din  -  2048 

The  relationship  between  Vin  and  In  is; 


In=204.8  X  Vin 


(72) 


In  our  digital  calculations,  Yi  is  the  measurement 
represented  by  In.  since  Yi;  -100<Yi<+100  ,  the  range  of  In 


IS 


In:  -100<ln<+100 


This  restricts  Vin  to; 


Vin:  -.4883v  <  Vin  <  +.4883v 


(73) 


Note  that  Yi=100  represents  some  scaling  factor(SF)  times  the 
position  in  nondimens ional  length  units  (NL=61.5").  The 
maximum  measured  deflection  is  |.0615“|  yielding: 


.0615"  or  .001  NL 
SF  =  lOexpS 


So  using  the  maximum  voltage  of  .4883v  we  obtain  a  desired 
sensitivity  of: 


.4883/. 001  (V/NL)  =  488.3  (V/NL) 
488.3(V/NL)  X  1000  (mv/V)  =  4.883  x  lOexpS  (mv/NL) 
4.883  X  lOexpS  (mv/NL)/61.5  (in/NL) 

7939.84  peakmv/in 


(75) 


The  accelerometer  sensitivity  is  calibrated  asi 


11.69  mv/g  at  lOhz 


11.99  mv/g  at  30hz 


Since  the  first  and  second  natural  frequencies  have  the 
same  order  of  magnitude  contributions  at  the  chosen  sensor 
location,  a  sensitivity  of  11.84  peak  millivolts  per  g  will  be 
assumed.  The  acceleration  signal  is  fed  to  the  vibration 
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meter.  Since  the  meter  requires  an  input  sensitivity  of  10 
mv(rms)  per  g(rms),  a  gain  of  1.194  is  applied  to  the 
acceleration  signal  before  the  integration  takes  place.  After 
proper  calibration,  the  meters  display  acceleration,  velocity, 
and  displacement,  and  the  peak  to  peak  displacement  signal  is 
displayed  on  an  rms  voltmeter.  The  desired  sensitivity  of 
5614  mv(rms)  per  inch  (peak  displacement)  was  obtained  by 
attenuating  the  displacement  signal  by  .1671. 

The  output  of  the  digital  to  analog  board  is  proportional 
to  the  integer  value  sent  to  it  (Vout).  The  linear 
relationship  isthe  same  as  for  the  analog  input.  We  recall 
that  since  the  scale  factor  is  lOexpS  and  therefore  the  output 
is  internally  divided  by  128  to  keep  within  the  allowable 
analog  output  integer  range  OUT;-2048<OUT<+2048  ,  the  force 
desired  is  in  Ibf  x  781.25.  Therefore,  if  the  value  (OUT) 
computed  is  78, the  desired  force  is  0.1  Ibf.  A  desired  force 
of  .00128  Ibf  corresponds  directly  to  10/2048  volts,  or 
4.8828  mv.  The  D/A  board  output  sensitivity  is  therefore: 

4.8828  mv 

-  —  or  3814.7  mv/lbf  (77) 

.00128  Ibf 

The  impedance  head  between  the  shaker  and  beam  measures 
force  with  a  sensitivity  of  557  mv(peak)/lbf .  The  following 
diagram  functionally  depicts  the  sensitivity  matching  which 
must  occur  prior  to  force  input  to  the  beam  after  being 
attenuated  by  a  factor  of  PS3  (potentiometer  #3  setting). 


3814.7mv/lbf — ►  PS3 


must  be 
557  mv/lbf  + 


Shaker  and 
Beam 
Dynamics 


557  mv/lbt 


Force 

►Applied 


Figure  3.  Actuator  Sensitivity  Matching. 


From  Figure  3  it  is  clear  that  the  output  force  signal 


must  be  attenuated  by: 


557 

PS3  =  -  =  0.1460 

3814.7 
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RESULTS/CONCLUSIONS 


Modal  Identification.  The  results  of  the  sine  dwell 
modal  identification  (Figures  4  through  7)  and  the  spectrum 
analysis  (Figures  8  through  10)  are  compared  to  the  computed 
values  in  Table  I.  The  ability  to  consistantly  excite  the 
beam's  first  mode,  using  the  sine  dwell  technique,  was  limited 
by  the  frequency  drift  in  the  function  generator.  Also, 
nonlinearities  were  discovered  when  too  much  force  was  applied 
by  the  shaker.  As  expected,  when  over  excited,  the  response 
peak  is  nolonger  symmetric.  Also,  Figure  11  graphically 
depicts  how  approaching  the  natural  frequency  from  lower 
frequencies  yields  drastically  different  results  than  when 
approached  from  above  when  nonlinearities  are  present. 

There  are  basicly  four  modeling  errors:  the  unmodelled 
accelerometer  mass  (2gm),  shaker  core  mass  (341gm),  the 
inability  to  obtain  perfectly  fixed  end  conditions,  and  the 
passive  damping  due  to  the  shaker.  The  first  three  will  lower 
the  measured  natural  frequencies  as  shown  in  Table  I.  The 
damping  due  to  the  shaker  supports,  however  was  not  noticeable 
until  the  control  experiment  results  were  examined. 


10.5 

11.0  11.5  Hz 

30 

31  32  33  34  35  Hz 

Figure 

4 .  Mode  1 

Figure  5.  Mode  2 

Sine  Dwell 

laentif ication. 

Sine 

Dwell  Identification. 

65  66  67  68  69  Hz 


111  115  119  Hz 


B'igure  6.  Mode  3 
Sine  Dwell  Identification. 


Figure  7.  Mode  4 
Sine  Dwell  Identification 
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Table  I 


Modal  Parameters. 


Natural 

Frequencies-to  (Hz) 

1 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Computed 

13.63 

37.57 

73.65 

121.74 

Sine  Dwell 

11.2 

32.55 

66.7 

115.1 

Spectral 

Analysis 

11.2 

35.3 

69 

109 

Assumed  in 

Model 

13.63 

37.57 

73.65 

121.74 

Damping  Ratio  (X  ) 

.L 

Sine  Dwell 

.0053 

.0031 

.0019 

.0043 

Assumed  in 

Model 

.005 

.005 

.005 

.005 

Software  Development.  The  main  emphasis  in  the  software 
developed  was  decreasing  the  cycle  time  tor  each  control 
output.  Both  the  measurement  sampling  and  the  control  output 
were  performed  in  the  same  program  loop.  Thus,  sampling  rate 
drove  both  the  observer  and  controller  sensitivity.  For  a  two 
mode  model,  it  was  desired  to  obtain  at  least  5  samples  per 
cycle.  Signed  16  hit  integer  arithmetic  kept  the  computation 
time  down.  The  system  was  upgraded  from  2.5mHz  to  5mHz 
operation  and  assembly  language  instructions  were  optimized. 
Only  four  samples  per  cycle  on  the  second  mode  were  achieved 
at  a  sampling  rate  of  125Hz  on  the  optimized  assembly  language 
subroutine  (Appendix  5).  When  it  was  determined  that  the 
sampling  rate  could  not  be  increased  further  without  obtaining 
another  processor  or  going  to  a  single  mode  approximation,  the 
two  mode  control  was  attempted. 


The  data  obtained  from  the 


static  deflection  testing  was  stored  on  disk,  then  plotted  (as 
in  Figures  19  through  22)  and  tabulated  in  Appendix  5.  The 
spectral  analyzer  output  a  plot  of  the  acceleration  PSD.  The 
half  power  method  was  used  again  to  compute  the  measured 
damping  for  each  mode.  The  expected  pole  migration,  obtained 
from  the  eigenvalues  of  Eq  42,  and  the  pole  migration  actually 
obtained  is  presented  in  Figures  12  and  13.  The  differences  in 
measured  natural  frequencies  seen  in  Figures  14  through  18  as 
compared  to  those  calculated  in  the  modal  parameter  section 
are  due  to  the  additional  shaker  used  to  input  the  noise  the 
noise  disturbance.  The  passive  damping  by  the  shaker  is 
evident  in  the  root  migration  tor  each  mode.  The  computed 
movement,  including  the  effect  of  the  residuals,  produced  only 
slightly  stable  roots,  while  the  actual  roots  became  more 
highly  damped  with  the  force  loop  closed  around  the  shaker. 
Note  that  no  change  occured  for  the  first  mode  between  the 
force  loop  closed  condition  and  the  control  applied  condition. 
The  second  mode  became  more  lightly  damped  however.  The 
residual  modes  were  expected  to  cause  a  more  lightly  damped 
set  of  modes  of  vibration,  but  their  actual  effect  is  masked 
by  this  passive  damping.  The  third  mode  damping  increased 
slightly,  and  the  fourth  mode  was  significantly  damped.  The 
damping  of  the  third  and  fourth  modes  was  not  due  to  the 
control  applied,  and  Figures  14,  15  and  16  clearly  demonstrate 
that  it  was  due  to  the  shaker  dynamics. 


Figure  12.  System  Pole  Movement 
Modes  1  and  2. 
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• 
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A 
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Figure  13,  System  Pole  Movement 

Modes  3  and  4. 


.2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  sec 

Figure  19.  Mode  1  Estimate  Force  Loop  Closed. 


•  2  .4  .6  .8  1.0  1.2  1.4  1.6  1.8  sec 

Figure  20.  Mode  1  Estimate  Control  Applied. 
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The  damping  of  the  third  and  fourth  modes  was 
anticipated  as  Herrick  (Ref  3)  discovered  similar  problems 
where  the  shaker  acted  as  a  passive  dashpot  when  used  in  this 
capacity.  He  found  this  to  be  the  result  of  the  shaker's 
flexure  supports  and  restored  the  beam's  resonant  behavior  by 
a  positive  velocity  feedback  loop  around  the  shaker. 

The  time  responses  of  the  estimate  of  mode  1  (Fig  19  and 
20)  for  0=200  show  that  the  more  lightly  damped  higher  modes 
predominate  shortly  after  control  is  applied.  Similarly,  for 
the  estimate  of  mode  2  (Fig  21  and  22),  the  second  mode 
quickly  predominates  and  a  much  cleaner  estimate  of  the  modal 
amplitudes  is  produced  for  this  reason. 

In  summary,  the  effect  of  this  controller  was  to 
slightly  damp  the  first  mode  and  negatively  damp  the  second. 
The  third  and  fourth  modes  were  damped,  however  this  is 
attributed  to  the  passive  damping  of  the  shaker.  Limitations 
on  the  control  applied  include  sampling  rate,  truncation 
errors  due  to  scaling  and  use  of  integer  arithmetic,  a  lack  of 
a  positive  velociy  feedback  loop  around  the  shaker,  and  a  lack 
of  phase  compensation. 
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RECOMMENDATIONS 


ff* 


The  experience  gained  in  implementation  of  optimal 
control  algorithms  in  this  experiment  has  helped  to  identify 
some  sensitive  areas  which  may  affect  performance. 

The  sensor  and  actuator  locations  are  very  important. 
Although  the  sensor  did  well  at  the  location  chosen,  the 
actuator  did  poorly  due  to  its  lack  of  control  authority  over 
the  first  mode.  The  optimum  actuator  position  is  dependent  on 
the  relative  cost  associated  with  not  controlling  each  mode. 
Therefore,  for  this  experiment  the  actuator  should  be  located 
nearer  the  center  of  the  beam.  In  addition,  the  diagonal 
elements  of  [Q]  should  be  weighted  to  reflect  these  costs  when 
calculating  the  optimum  controller  and  observer  gain 
matricies. 

Shaker  damping  was  a  significant  factor  inthe  degraded 
controller  performance  and  should  be  compensated  for  using  a 
velocity  feedback  path  around  the  shaker.  An  impedance  head 
which  outputs  both  acceleration  and  force  simultaneously  is 
ideal  for  this  implementation. 

The  need  for  faster  computation  speed  is  obvious;  however 
no  matter  what  hardware  constraints  are  imposed,  the  control 
algorithm  will  have  to  be  optimized  to  yield  acceptable 
computation  rates  and  accuracy.  A  working  knowledge  of  the 
specific  machine  code  utilized  will  be  required  for  this. 

Finally,  new  methods  for  exciting  the  beam  without  the 
addition  of  a  second  actuator,  such  as  acoustic  vibration  and 
the  use  of  a  load  cell  equipped  hammer,  should  be  examined. 
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Appendix  1 


^  □  * 


^‘D  o  □ 


o  □ 


m 


shaker.. . Detune  Clip-Full  Up. 

Low  Frequency  Function  Generator ..  .Function-Sind,  Ainp<10. 

Power  Amplifier . Cathode-1. 

DC  Power  Supply . H.V.-ON,  DC  Amps=200. 

True  RMS  Voltmeters . Use,  Meter. 

Universal  Counter-Timer . . . AC  Input,  Func. -Period 

Time-  N>10exp2. 

Accelerometer. . . .  .&Blect  .Iso.  Mechanism. 

battery  Power  Unit.... . Gain^lOO. 

Impedance  Head..... . Use  Force  Output. 

Electrical  Connection . Beam  Center-Modes  1&3. 

Offset-Modes  2&4. 


Figure  23.  Sine  Dwell 
Modal  Parameter  Identification 
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0  O  □ 


□  on 


f  F  •••••••••••••••••••••••••••  aSfiinG  ds  foir  PiQuiTG  23  • 

L  -  Noise  Generator . Seq.Leny th-ln£initer  Gaussian  Output. 

M  -  Correlator . .Punc.-Autocor, , Averaging-Exponential , 

Dither-ON# Delay  Ottset-0, Interpolation-On,Input-A-AC. 

N  -  Spectrum  Display . Input-1 ,Disp-Spectrum, Real ,Window-of£ 

Gain-5 ,  interpolation-Autd, Marker-off  »Vertical-Moclulus 
Horizontal-preq, , Phase  shift-0 ,Cal-Off , Record-scan, 

0  -  X-Y  Recorders. . .Time  Base-Otfr Sweep  Rate-5 , Line-On, 

Pen-1 i ft, Servo-On ; Polar i ty-X : +RT , Y :+UP ; Reponse-X : S low 
YjSlow;Range-X:lv/in,Y: ,5v/in. 


Figure  24.  Spectrum  Analysis 
Modal  Parameter  Identification 
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Figure  25  Legend. 


A  -  Video  Terminal. 

B  -  8"  Floppy  Disk  Drive. 

C  -  S-100  Mainframe;  1  -  SBC-100;  Serial  Input/Output  to  A. 

2  -  Disk  Drive  Board. 

3  -  Memory  Board. 

4  -  A/D  Board:  Input  on  Chan  0  from  Ql . 

5  -  D/A  Board;  Chan  0  -  Q3 . 

Chan  1  -  X  on  plotter. 

Chan  2  -  Y  on  plotter. 

Chan  3  -  PenUp  on  plotter. 

D(l&2)  -  Shaker . .....Same  as  Figure  23. 

E-  PortableCabinet;SameasFigure23, except  E2-Input  from  Q2 

Out  Disp.->01 
Connect  to  03 

HI  -  Accelerometer. 

H3  -  Battery  Power  Unit . Gain=10. 

I  -  Impedance  Head. 

L  -  Noise  Generator . Same  as  Figure  24. 

M  -  Correlator . Same  "  "  ". 

N  -  Spectrum  Display . . . Same  "  "  ". 

O  -  X-Y  Recorder;  Continuous  Vibration  Input  -  Input  from  N. 

Initial  Displacement  Input  -  Input  from  C5. 

0  -  Analog/Hybrid  Computer;  Pot.  Settings;  01=  .1671 

02=  .1194 
03=  .1460 
04=  .2571 

R  -  Initial  Condition  Displacement  Mechanism;  Rod  w/wall  pivot 
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LABORATORY  EQUIPMENT 


a 

1  o 


I 


« 


i 


A.  Heathkit  Terminal(#5) :Model  H19-A  -  SN  17-43722. 

B.  Tarbell  Floppy  Disk  Drive;Model  VDS-TTD  -  SN  2267/0482. 

C.  ECT( Electronic  Control  Technology ) Mainframe  MB  -10. 

1.  SBC-100  Single  Board  Computer  -  SD  Systems 

2.  Single  Density  Floppy  Disk  Interface  Board 

Tarbell  -  SN  10867. 

3.  Expanaoram  l i i-Expandable  Random  Access  Memory- 

SD  Systems  -  SN  AW  0100160. 

4.  AIM-12,  12  bit  Analog  Input  Module-Dual  Systems 
Corporation. 

5.  AOM-12,12  bit  4  Channel  S-100  Digital  to  Analog 
Converter  -  Dual  Systems  Corporation  -  SN  1039. 

D.  Shakers  -  MB  Electronics  Vibration  Test  Equipment 

1.  Model  SAC  -  SN  186. 

2.  Model  SDA  -  SN  271. 

E.  Portable  Cabinet 

1.  Accelerometer  Integrator/Amplifier 
MBElectronics-Model  N504-SN  PS4110  and  PS4111. 

2.  Vibration  Meter  -  MB  Electronics  -  Model  N499 

SN  168  and  SN  151. 

3.  Vibration  Indicator-  MB  Electronics-  Model  N497 
SN  162  and  SN  173. 

4.  Low  Freq  Function  Generator  -  Hewlett  Packard 
Model  202A  -  SN  325-13268  and  4156. 

5.  125VA  Power  Amplifier  -  MB  Electronics  Model 
2120MB  -  SN  904  and  913. 
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6.  DC  Power  Supply  -  Hewlett  Packard  -  Model712B 
SN  002-05035  and  SN  002-05037. 

SN  002-05082  and  SN  002-05083. 

True  RMS  Meters  -  Ballantine  Laboratories  -  Model  320 

1.  SN  1537(EE  Lab). 

2.  SN  1538(EE  Lab). 

Universal  Counter-Timer  -  Computer  Measurements  Co. 
Model  726C  -  SN  49300/4H5581. 

I.C.P.  Accelerometer-PCB  Piezotronics  Inc. 

1.  Model  303A  -  SN  1935. 

2.  Model  303A  -  SN  1934. 

3.  Battery  Power  Unit  -  PCB  -Model  480C06  -SN  1247. 

4.  Battery  Power  Unit  -  PCB  -Model  480D06  -SN  2062. 
Impedance  Head  -  Wilcoxon  Research 

1.  Model  Zll  -  SN  900. 

2.  Model  Zll  -  SN  802. 

Oscilloscope  -  Tektronics  Inc.-  Model  116E-  SN  002800. 
DC  Power  Supply  -  H.P.  -  Model  6205B  -  SN  1949A16663. 
Noise  Generator  -  H.P.  -  Model  3722A  -  SN  1722U01822. 
Correlator  -  H.P.  -  Model3721A-  SN  1112A00104. 
Spectrum  Display  -  H.P.  -Model  3720A  -  SN  1215U00120. 
X-Y  Recorders  -  H.P.  -  Model  7045B  -  SN  2129A00481. 
S-100  Bus  Extender  Board-Bob  Mullen  Computer  Products 
Model  TB4  -  SN  0281. 

TR-48  Analog/Hybrid  Computer-EAI-  Model  45.117-SN  568. 
Initial  Condition  Displacement  Mechanism;  Rod. 


<D 


Where 

V  (x)  =  C  sin(y8  x)  -  C  sinh(^  x)+Z  C  [cos(yQ  x)-cosh{^  x)]  (14) 

i  li  1  iil  i  i 


For  L  =  1,  and  A  =  4.73  =>  Z  =  -1.017809 

1  1 

=>  C  =  .80452 
1 

=>  C  *  -.81885 
2 


(h  =  7.853  =>  Z  =  -.99922 

2  2 

=>  C  =  .819261 
1 

=>  C  *  -.818622 
2  ' 


/3  =  10.996  =>  Z  =  -1.00003 
3  3 

=>  C  =  .818953 
1 

=>  C  =  -.81898 
2 


^  *  14.1372=>  Z  =  -.999998 
4  4 

=>  C  =  .827169 
1 

=>  C  =  -.827168 
2 
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So  that  the  normalized  eigenvectors  are; 


V  (x)  =  .80452  [sin(4 .73x)  -  sinh{4,73x)] 

1 

+  .81885 [cosh(4.73x)  -  cos(4.73x)]  (17) 

V  (x)  =  .81926 [sin(7.853x)  -  sinh( 7 . 853x ) ] 

2 

+  .81862 (cosh(7.853x)  -  cos(7.853x)]  (18) 

V  (x)  =  .81895 [sin(10.996x)  -  sinh( 10 . 996x ) ] 

3 

+  .81898 [cosh(10.996x)  -  cos( 10 . 996x ) ]  (13c) 

V  (x)  =  .82717 [sin(14.1372x)  -  sinh( 14 . 1372x ) 1 

4 

+  .82717 (cosh(14.1372x)  -  cos( 14 . 1372x ) ]  (13d) 


The  values  of  these  mode  shapes  are  then: 


Sensor(x  =  .852698  NL)  Actuator(x  =  .079956  NL) 


s 

a 

V  (X) 

1 

.3071239 

.1026195 

V  (x) 

2 

-.6945324 

.2539699 

V  (x) 

3 

1.0171641 

.448457 

V  (x) 

4 

-1.229404 

.6626658 

ADoendix 


This  appendix  contains  tor  Q  =  20,  200,  500,  6000: 

A.  OPTCON  generated  controller  and  observer  gains  [G] 

T 

and  [K]  (labeled  feedback  matrix  in  computer 
output) . 

B.  Closed  loop  system  matrices  [ [Ac]  +  [Be]  [G] - [k]  [Cc]  ] 
without  residuals  present  ana  their  eigenvalues. 

C.  Eigenvalues  of  the  closed  loop  system  matrices  (with 

residuals  present)  of  the  form: 


[  [AcJ  +  [Bc]  [G]]  [BC]  [G]  [0] 

[0]  [ (Ac]-[K] [Cc] ]  [K] [Ccr] 

[Bcr] [G]  [Bcr] [G]  [Acr] 


D.  Computed  [H] ,  [F],  and  [G]  matrices  for  computer 
implementationof  TL  =  [H]^+  [F]Y  and 


STEADY  STATE  FEEDBACK  GAIN  FOF<  THE  LINEAR  REGULATOR 
DX/DT  =  A#X  +  B*U 
2J  =  INTEGRAL  OF  X'i*:«*X  f  U'*R*U 

WHERE 


A 


0* 

0. 

l*0000E+00 

0. 

0* 

0. 

0. 

1 .OOOOEiOO 

7.3321E+03 

0. 

--8*5630E-"01 

0. 

0* 

-5.57()9EI04 

0» 

-2.3603E400 

D 

0. 

0. 

1*0262E-01 

2.5397E-01 


0 


2.0000E401 

0» 

0* 

0. 

0« 

2.0000E+01 

0. 

0. 

0, 

0, 

2.0000E{01 

0. 

0. 

0» 

0. 

2*0000E401 

R 

l.OOOOEfOO 


RICCATI  SOLUTION 

8.0238E+04  -7.1770E-01  1.3A37E-03  -1.7360E-01 
7.1770E-01  2.2375Ei05  1.3.190E+00  1.7934E-04 
1.3637E-03  1.3190E+00  1.0943ET01  -3.5887E-05 
1.7360E-01  1.7934E04  -3.5887E-05  4.0164E+00 


FEEDBACK  MATRIX 

-4.3949E-02  1.3540E-01  1.1230Ei00  1.0200E+00 


CLOSED  LOOP  EIOENUALUES 
REAL  IMAGINARY 

-  .48577E+00  •-  .85626E+02 
-.48577E+00  ♦85A26Ef02 
-.13097E+01  -.23602E+03 
-.13097E+0I  ♦23602Ef03 
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STEADY  state:  FEiEDBACK  GAIN  FOR  THE  LINEAR  REGULATOR 


BX/DT  =  A*X  +  B*U 
2J  =  INTEGRAL  OF  X'*a*X  +  Lr»R*U 

WHERE 

A 


0. 

0. 

-7,3320E+03 

0. 

0. 

0. 

-5.5709E+04 

1 .OOOOE+00 

0. 

-8,5630E“01 

0, 

0. 

1 .OOOOE+OO 

0. 

-2.3603E+00 

B 

♦0712E-01 

.9453E-01 


Q 


OOOOE+01  0, 

0. 

0* 

2 ,0000  fOl 

0» 

0. 

0« 

2»OOOOE+Ol 

0* 

0. 

0. 

2»OOOOE+Ol 

R 


.OOOOE+00 


RICCATI  SOLUTION 

.0830Ef00  2»7947E-04  -6,9191E+00  -6»3415E+00 
7947E-04  3.1945E+00  8.3428E--01  ~7.5389E+00 
9191E+00  8.3428E-01  5»9252E+04  8.7636E+00 
3415E+00  -7,5389E+00  8*7636Ei00  1 . 7794E+05 


FEEDBACK  MATRIX 

2.4823E+00  -2.2186E+00  -2,7045E+00  3.2884E{00 


CLOSED-LOOP  EIGENVALUES 
REAL  IMAGINARY 


- .80934E+00 
-.80934E+00 
- ♦19506E+01 
-.19506Ef01 


-»85623E+02 

.85623E+02 

-.23602ET03 

.23602E+03 
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Y 


EAI: 


STATE  rEEDBACK  GAIN  EDFi;  THE  LINEAR  REGULATOI 


DX/DT  -  A*X  {  i:{>iiU 


2J  ■•=•  INTEGRAL  OF  X'XcOilcX  +  U')f:R*U 

WHERE 


C.  0. 

0 .  0 . 

/■*3321Li03  0. 

0*  ■■5.t;>709Ej04 


l*000OL{00  0. 

0.  :l,  .OOOOEiOO 

B*5630E-  C;l.  0. 

0*  -2.3&03E{00 


B 


0. 

0* 

;U02u2E  01 
2*S397E-0:l 


Q 


nr 


2.0000Ef02 

0, 

0  t. 

0. 


0. 

2.0000Li02 

0  * 

0» 


0  ♦  0 , 

0  *  0 « 

2vOOOOEi02  0» 

0.  2.0000LI02 


R 


1 »OOOOEiOO 


S,  7709!:}  OS 
■A.43:L3E}0;1. 
:l.  .32S6E-  02 
■  V  >  3 A  2 1.;.  {  0  (/ 


r  LED  BAG  I- 
•  2.3<?.90Lv00 


DI...UTION 

6.4313E-I01 
1.  .673SE  {06 
7.0939EI01 
i  ♦3I94E-  03 


MATRIX 
7. 200  IE { 00 


1 .325AE-02 
7 *09390101 
7.8710ET01 
-3.3003E- 03 


0.0764E : 00 


■••9.33(!>2Ei00 
1  *3194E--03 
■-3.3003E-03 
3.0039i;:{01 


7.32CnE{00 


CLOSED-  E.OnR 
REAL. 

♦  21409L-IC1 
*21-^:0VL  !01 
.OISSALIOO 
-.S42SAM  00 


EIGEF’OAIUES 
I MAO 1 NARY 

.23302E! 03 
.23602E{ 03 
.05624EI02 
.0SA24EI 02 
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STEADY  STATE  FEEDBACK  GAIN  FOR  THE  LINEAR  REGULATOR 


DX/DT  =  A*X  +  B*U 
2J  =  INTEGRAL  OF  X'*Q*X  +  U-'JtcRJlcU 

WHERE 

A 


0. 

-7.3321E+03 

0. 

0. 

0* 

-5 « 

5709E+04 

OOOOE+00  0, 

-8.5630E-01 

0* 

1 .OOOOE+00 

0. 

-2. 

3603E+00 

B 

.0712E-01 

.9453E-01 


Q 


2.0000E+02 

0, 

0, 

0* 

0, 

2.0000E+02 

0. 

0* 

0, 

0, 

2»0000E+02 

0. 

0. 

0. 

0. 

2.0000E+02 

R 

l»OOOOEfOO 


R ICC AT I  SOLUTION 

3,7890E+01  1,7737E“02  -3,2434E+01  -1.4915E+02 

l,7737E-02  1.6054Ef01  1,9598E+01  -3.7899E+01 

-3.2434Ei01  1.9598E+01  2»7751E+05  4.0954E+0 

-1.4915E+02  ~3.7899E+01  4.0954E+02  8.9448E+0 


FEEDBACK  MATRIX 

1.1625E+01  -1.1145E+01  ~2.3573E+01  -1.94P5E+01 


CLOSED- LOOP  EIGENVALUES 
REAL  IMAGINARY 

-.22134E+01  -.85599E+02 
-.22134E+01  .85599E+02 
-♦50501E+01  “.23597E+03 
-.50501E+01  .23597E+03 
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STEADY  STATE  FEEDBACK  GAIN  FOR  THE  LINEAR  REGULATOR 


DX/DT  =  A*X  {  B*U 
2J  =  INTEGRAL  OF  X'*Q*X  i  U'*R*U 

WHERE 

A 

♦  0.  l»0000E+00  0. 

.  0.  0.  l.OOOOE+00 

♦3321E+03  0.  -8»5630E-01  0. 

»  “5.3709ET04  0.  -2»3603E-I00 


B 


!o2&2E-01 

.5397E-01 


Q 


5.0000E+02 

0. 

0. 

0. 

0. 

5.0000E+02 

0. 

0, 

0. 

0. 

S.OOOOE+02 

0. 

0, 

0. 

0. 

'  5.0000E+02 

R 

l.OOOOE+00 


RICCATI  SOLUTION 

1.1092E+06  -3,4741E+02 
3,4741E+02  3.2737EI06 
2,8635E-02  2*6663Ei02 
3,5088E+01  -2,2315E-03 


FEEDBACK  MATRIX 

8.9084E+00  2.7361EJ01  1.5522E+01  1.4921E+01 


CLOSED-LOOP  EIGENVALUES 
REAL  IMAGINARY 

-♦12246E+01  --,85619E+02 
“.12246E+01  .85619E+02 
-♦30748Ef01  -♦23601E+03 
-.30748E+01  »23A01E+03 


2.8685E-02  -3t5088E+01 
2.6663C;f02  -2.23151  03 
1.5130Ef02  -1.795C:  02 

-1.7958E-02  5.8757E+01 
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STEADY  STATE  FEEDBACK  GAIN  FOR  THE  LINEAR  REGULATOR 


DX/DT  =  A*X  f  B*U 
2J  =  INTEGRAL  OF  X'*0*X  +  U'*R*U 


WHERE 


0. 

A 

0. 

~7.3321E+03 

0. 

0. 

0. 

0. 

-5.5709E+04 

l.OOOOE+00 

0. 

-8.5630E-01 

0. 

0. 

l.OOOOE+00 

0. 

-2.3603E+00 

D 

3.0712E-01 

-6.9453E-01 

0. 

0. 


Q 


5.0000E+02 

0. 

0. 

0. 

0. 

5.0000E+02 

0. 

0. 

0. 

0. 

5,0000E+02 

0. 

0. 

0. 

0. 

5.0000E+02 

R 

l.OOOOE+00 


RICCATI  SOLUTION 

6,4446E+01  8.0944E-02 

8.0944E-02  2»7694EI01 

-5»5234E+01  5.7243E-f01 

-4»3637E+02  -6 . 5498E+01 


FEEDBACK  MATRIX 

1.9737E+01  -  1.9209E+01  ■-5»6720E+01  -a.8523E+01 


CLOSED- LOOP  EIGENVALUES 
REAL  IMAGINARY 

-.34594E+01  -,85558E+02 
-.34594E+01  ,85558E+02 
-.78S04E+01  -.23590E+03 
-.78504E+01  ,23590E+03 


-5.5234E+01  -4.363'E+02 
5.7243E+01  -6.549  +01 

4.7136E+05  1.73:  1+03 

1.7321E+C3  1.544  CA06 


STEADY  STATE  FEEDBACK  GAIN  FOP;  THE  LINEAR  REGULATOR 


DX/DT  =  A*X  +  B*U 
2J  -  INTEGRAL  OF  X'»03l:X  +  U'*R*U 

WHERE 


A 


0. 

0. 

l.OOOOE+00 

0. 

0. 

0. 

0. 

l.OOOOE+00 

7.3321EF03 

0. 

-8.5630E~01 

0. 

0. 

••-5.5709E+04 

0. 

-2.3603E+00 

B 

0. 

0. 

1 .0262E-01 
2.5397E-01 


0 


6,0000Ef03 

0. 

0. 

- 

0. 

A.OOOOE+03 

0. 

0. 

0. 

^>.0000EF03 

0. 

0. 

0. 

6*0000E+03 

1 .OOOOEFOO 


RICCATI  SOLUTION 

4.9752E+06  --2.3089Et04  -1.8803E+00  -7*2070Ei02 
-2.3089E+04  1.5129E+07  5.4899E+03  -2.7877E+00 
■•1.8803E+00  5.4899E  +  03  6.8027E+02  -1.2009E  +  00 
-7,2070E+02  -2,7877E{00  -1.2009E+00  2.708?  f02 


FEEDBACK  MATRIX 

-1.8323Ef02  5.6267Ei02‘  6.9567Ef01  <!>,,8671E+0l 


CLOSED- LOOP  EIGENUALUES 
REAL  IMAGINARY 

•  .  399r.4E  f01  -♦85535E+02 

.399S4E+01  .85535E+02 

~.9B995E+01  -■  .23582E+03  , 

-.98995E+01  .23582E+03 


[  [AC]  +  [BC]  [G]  -  [K]  [Cc]  ] 


For  0=20:  Eigenvalues: 


-.7624 

1.724 

1.000 

0 

-.7517  +/- 

J 

85.62 

.6814 

-1.541 

0 

1.000 

-7331. 

-1.864 

-.7411 

.1047 

-1.821  +/- 

J 

236.0 

-1.021 

-55710 

.2852 

-2.101 

For  0= 

200: 

-3.572 

8.077 

1.000 

0 

-1.800 

+/- 

J 

85.63 

3.424 

-7.744 

0 

1.000 

-7325. 

-15.62 

-.0275 

.7829 

-4.083 

+/- 

J 

235.9 

5.384 

-55720 

2.051 

-.4228 

For  0= 

500: 

-6.062 

13.71 

1.000  0 

-2.666 

+/- 

J 

85.65 

5.899 

-13.34 

0  1.000 

-7316. 

-36.58 

.7366  1.531 

-5.952 

+/- 

J 

235.8 

24.93 

-55760 

3.942  1.429 

For  0= 

6000: 

-22.83 

51.64 

1.000 

0 

-8.427  +/-  J 

86.42 

22.62 

-51.16 

0 

1.000 

-7157. 

-379.7 

6.277 

7.047 

-17.89  +/-  J 

232.1 

484.6 

-56770 

17.65 

15.08 

The  eigenvalues  of  the  closed  loop  system  with  resid 
present  are; 


E'or  0=20: 


-.8096 

+/- 

J 

85. 

-.3706 

+/- 

J 

85. 

-1.9504 

+/- 

J 

236 

-1.0505 

+/- 

J 

236 

-2.3139 

+/- 

J 

462 

-2.3139 

+/- 

J 

764 

-0.0137 

+/- 

J 

85. 

-2.2145 

+/- 

J 

85. 

-0.2115 

+/- 

J 

236 

-5.0518 

+/- 

J 

236 

-2.3135 

+/- 

J 

462 

-3.8249 

+/- 

J 

764 

For  0=500; 


+0.3684 

+/- 

J 

85. 

-3.4605 

+/- 

J 

85. 

+0.7140 

+/- 

J 

235 

-7.8489 

!/- 

J 

235 

-2.3114 

+/- 

J 

462 

-3.8266 

+/- 

J 

764 

For  0=6000; 


+03.1510 

+/- 

J 

85 

-11.9052 

+/- 

J 

84 

+07.4927 

+/- 

J 

235 

-26.7056 

+/- 

J 

233 

-02.1238 

V- 

J 

463 

-03.9725 

+/- 

J 

764 

I 


I 


6224 

6263 

.0198 

.0251 

.7612 

.9125 


6291 

6019 

.0505 

.9636 

.7740 

.9085 


6397 

5716 

.9869 

.8623 

.8053 

.8986 


.8333 

.9542 

.1097 

.8435 

.5620 

.6554 


0=20 


[H] 


•7697  .00572 

.002269  -.3112 

-53.85  -.4197 

-.8969  -224.8 


007345 

00002155 

7699 

.001568 


00004066 

004035 

0005213 

3135 


IF]  =  [  .01812  -.008787  -.5747 
[G]  =  [-.043949  .1354  1.123 

Pseudostate  Form; 


2.900  ] 

1.020  ] 


[H]  ' 


.7697 

.002269 

-.5385 

-.008969 


.00572 

-.3112 

-.004197 

-2.248 


.7345 

.002155 

.7699 

-.001568 


.004066 

.4035 

-.0005213 

-.3135 


[F] 

[G] 


I  .1812  -.08787  -.05747 

[-.00439  .01354  11.23 


Multiply  [H] '  and  [F] 


by  128  and  form  integer  matrices; 


[H]  = 

99 

1 

94 

1 

0 

-  40 

0 

52 

-69 

-  1 

99 

0 

-  1 

-288 

0 

-40 

IF]  =  [ 

23 

-  11 

-7 

3  ] 

[G]  =  [ 

-  0 

0 

11 

10  ] 

0=200 


pi  o 


[H]  = 


.7513 

.01078 

-53.38 

-4.493 


[F]  =  I  .08209 

[G]  =  [-2.3698 


.02427 

-.3220 

-2.316 

-220.9 


-.04310 

7.2801 


,007287 

,0001208 

,7765 

,005108 


•2.688 

8.0764 


.0001933 

.003964 

-.001432 

-.2936 


14.14  ] 

7.6288  ] 


Pseudostate  Form; 


[H] 


.7513 

.01078 

.5338 

.04493 


.02427 

-.3220 

-.02316 

-2.209 


[F] 

[G] 


[  .8209  -.431 

[-.23698  .7280 


,7284 

,01208 

,7765 

,005108 


-.2688 

80.764 


.01933 

.3964 

-.001432 

-.2936 


1.414  ) 

76.288  ] 


Multiply  [H]  and  [F]  by  128  and  form  integer  matricesi 


[HI  = 


[F] 

[G] 


96 

3 

93 

2 

= 

1 

-  41 

2 

51 

-68 

-  3 

99 

0 

-  6 

-283 

-  1 

-38 

T 

=  [105 

-55 

-34 

181] 

=  ( 

0 

1 

81 

76] 

0=500 


[H] 

.7358 

=  .01788 

-53.01 
-7.646 

.03832 

-.3320 

-4.159 

-217.9 

.007244 

.0002152 

.7834 

-.006418 

.0003274 

.003907 

-.001706 

-.2759 

IFJ 

=  [  .1353 

-.07268 

-4.551 

23.74 

[G] 

=  [-8.9084 

27.361 

15.522 

14.921 

Pseudostate 

Form; 

f 

I 

.7358 

.03832 

.7244 

.03274 

[H] 

.01788 

-.3320 

.02152 

.3907 

-.5301 

-.04159 

.7834 

-.001706 

-.07646 

-2.179 

-.006418 

-.2759 

9 

T 

IF] 

=  [ 

1.353 

-.7268 

-.4551 

2.374 

] 

[G] 

=  [ 

-.8908 

2.7361 

155.22 

149.21 

1 

Multiply  [H] 

and 

[F]  by 

128  and  form  integer 

matrices: 

9 

94 

5 

93 

4 

[H] 

2 

-42 

3 

50 

-68 

-  5 

100 

0 

-10 

-279 

-1 

-35 

Li» 

— 

$ 

[Fl 

[173 

-93 

-58 

304] 

[G] 

[ 

-  1 

3 

155 

1491 

Q=6000 


• 

6527 

.09193 

.007073 

.001147 

[H] 

= 

« 

05591  - 

.3879 

.0008353 

.003601 

-50.87 

- 

16.65 

.8357 

-.0007079 

-26.19 

•• 

203.9 

.007458 

-.1694 

T 

[F] 

= 

« 

4181 

— 

.2429  - 

16.69 

77.01 

[G] 

= 

-183.23 

562.67 

69.507 

68.671 

Pseudostate 

Form 

# 

.6527 

.09193 

.7073 

.1147 

1 

[H] 

= 

.05591 

-.3879 

.08353 

.3601 

-.5087 

-.1665 

.8357 

-.0007079 

-.2619 

-2.039 

.007458 

-.1694 

J 

IF]  ' 

= 

4.181 

-2.429 

-1.669 

7.701 

[G] 

as 

[-18.323 

56.267 

695.07 

686.71 

Multiply  [H] 

f  1 

and  [F] 

by 

128  and 

form  integer  matrices 

84 

12 

91  15 

IH]  = 

7 

-  50 

11  46 

-65 

-  21 

107  0 

-34 

-261 

1  -22 

T 

IF]  =  [535 

-311 

-214  986] 

With  Gain  =  100  at  accelerometer  amplifier: 

T 

IF]  =  15-3-2  10] 


Divide  by  18.323  here,  and  Multiply  by  18.323  outside 


icrete  Form  Derivation 


The  details  ot  the  discrete  torm  integrations  are 


provided  below: 


[A]  t  _  p  t  [A]  (t-^) 

X(t)  =  e  X(0)  +  \  e  [KlYC'l)  dT:  (50) 


Over  one  sample  period, T, 


_  [A]T  _  p  T  [A]  (T-T:) 

X(T)  =  e  X(0)  +  \  e  [K]Y(T)  d't  (50a) 


Assuming  Y('l)  is  constant  over  the  sample  period,  0  to  T, 


_  (A]T  __  lAjTpT  -[AR 

X(T)  =  e  X(0)  +  e  \  e  d^  [K]Y(0)  (50b) 


It  can  be  easily  shown  that: 


T  -[AlT:  -1  -lAj  T 

e  d't  =  -lA]  e 

0  t  =0 


(50c) 


Which  by  substitution  leaves: 


_  [AlT  _  lA]T  -1  -[A]T 

X(T)  =  e  X(0)  +  e  lA]  [[I]-e  ]tK]Y(0)  (50d) 


_  [A]T  _  -1  [A]T 

X(T)  =  e  X(0)  +  lA]  [e  -  [I]][K]Y(0)  (51) 


Appendix  5 


ontrol  Alaorithro  Computer  Listin 


The  computer  programs  used  for  control  implementation 
are  listed  below: 


1.  TASK9.  This  Fortran  80  main  body  program 
incorporates  the  storage  of  controller  matrices'  elements  and 
commands  the  subroutines  to: 

a)  Start  the  system  clock  counters ( STCLK ) . 

b)  Perform  the  control  function(CNTRL2M) . 

c)  Write  the  data  stored(time  response  of  the 
estimate  of  the  mode)  to  a  disk  f ile( WRTBUF ) . 

d)  Plot  the  discretized  time  response ( PLOT) , 
Support  subroutines  are: 

a)  Analog  to  Digital  and  Digital  to  Analog 
service  routines(AD  and  DA,  contained  in  ADIO). 

b)  A  simple  time  delay  routine  which  allows  the 
plotter  time  to  respond  to  directed  inputs(WAIT) . 


PROGRAM  TASK9 


PROGRAM  TASK9 
IMPLICIT  INTEGER(A-Z) 

C  DEFAULT  IS  INTEGER  *  2 

DIMENSION  F(4) ,G(4) ,HX(4) ,X(4) ,Y(iOOO) 

DIMENSION  IH(4,4) 

WRITE(1,22) 

22  FORMAT(4X, 'RUN  NUMBER( 1 , 2 , 3 ,4 ) ? ' ,/ ) 

READ( 1 ,24) RUN 

24  FORMAT(Il) 

C  THE  FOLLOWING  VALUES  COME  FROM  EQ  66  AND  68  IN  THESIS 
F(l)=105 
F(2)=-55 
F(3)=-34 
F(4)=181 
IH(1,1)=96 
IH(1,2)=3 
IH(1,3)=93 
IH(1 ,4)=2 
IH(2,1)=1 
IH(2,2)=-41 
IH(2,3)=2 
IH(2,4)=51 
IH(3,l)=-68 
IH(3,2)=-3 
IH(3,3)=99 
IH(3,4)=0 
IH(4,l)=-6 
IH(4,2)=-283 
IH(4,3)=-1 
IH(4,4)=-38 
G(1)=0 
G(2)=l 
G(3)=81 
G(4)=76 
DO  90  N=l,1000 

90  Y(N)=0 

DO  100  1=1,4 

X(I)=0 

HX(I)=0 

100  CONTINUE 

CALL  STCLK 

CALL  CNTL2M(F(1) ,G( 1 ) ,IH( 1,1 ) ,HX( 1 ) ,X( 1 ) ,Y( 1 ) ) 
CALL  WRTBUF(Y(1) ,1000, RUN) 

CALL  PLOT(Y(l) ) 

STOP 

END 


68 


SUBROUTINE  CNTL2M 


ENTRY  CNTL2M 
EXTERNAL  $AT,$M9 
CSEG 
CNTL2M: 

THIS  PROOGRAM  IS  DERIVED  FROM  THE  FORTRAN  COMPILER'S 
VERSION  OF  AN  EQUIVALENT  SUBROUTINE,  USING  TECHNIQUES 
AS  OUTLINED  IN  THE  THESIS  PAGE  8.  SOME  SUBROUTINES 
ARE  MADE  PROGRAM  SPECIFIC  TO  REDUCE  COMPUTATION  TIME. 
COMMENTS  ARE  DENOTED  BY  A  C  IN  THE  SECOND  COLUMN.  THE 
EQUIVALENT  S/R  IN  FORTRAN  IS  NOTED  BY  AN  F  IN  THE 
SECOND  COLUMN: 


SUBROUTINE  CNTL2M ( P ,G , IH ,HX ,X , Y ) 
IMPLICIT  INTEGER(A-Z) 

DIMENSION  Y(IOOO) ,HX(4) ,X(4) ,F(4) ,G(4) 
DIMENSION  IH(4,4) 

INITIALIZE  PROGRAM 


. RADIX 

16 

SHLD 

F 

XCHG 

SHLD 

G 

MVI 

A, 04 

LX  I 

H,IH 

CALL 

SAT 

;F 

DO  300  K=0,253 

;C  REMOVE  THE  NEXT  SEMI 

-COLON  FOR 

;  DOK : 

SHLD 

K 

;F 

DO  300  N=l,999 

,2 

LXI 

H,0001 

DON: 

SHLD 

N 

;F 

CALL  GETTIM(Y(N-H)  ) 

LHLD 

N 

DAD 

H 

XCHG 

LHLD 

Y 

DAD 

D 

SHLD 

TOO 

LHLD 

TOO 

.Z80 

IN 

A,(07AH) 

LD 

(HL) ,A 

IN 

A,  (07BH) 

INC 

HL 

LD 

(HL) ,A 

.8080 

;F 

OUT=0 

LXI 

H,0000 

SHLD 

OUT 

;F 

CALL  AD(VIN,0, 

80) 

LXI 

B, FIFTY 

LXI 

D,ZBRO 

LXI 

H,VIN 

LD 

(VALUE) ,HL 

LD 

( CHAN ) , DE 

LD 

(BASE) ,BC 

EX 

DE,HL 

LD 

A,(HL) 

LD 

HL,(BASE) 

LD 

C,(HL) 

OUT 

(C)  ,A 

INC 

C 

LD 

A,0 

OUT 

(C)  ,A 

DEC 

C 

IN 

A,(C) 

AND 

080H 

JR 

NZ,NRDY 

INC 

C 

INC 

C 

IN 

A,(C) 

LD 

B,A 

INC 

C 

IN 

A,(C) 

AND 

OFH 

LD 

D,A 

LD 

HL, (VALUE) 

LD 

(HL) ,E 

INC 

HL 

LD 

(HL) ,D 

2048 

LHLD 

VIN 

LXI 

D,0F800 

DAD 

D 

SHLD 

IN 

1=0,3 

LXI 

H,0000 

SHLD 

I 

J=0,3 

LXI 

H,0000 

SHLD  J 

230  HX(I+1)»IH(I+1,J+1)*X(J+1)+HX(I+1) 
LHLD  I 

DAD  H 

SHLD  TOO 

XCHG 

LHLD  HX 

DAD  D 

SHLD  TOl 

LHLD  J 

DAD  H 

DAD  H 

DAD  H 

XCHG 

LHLD  TOO 

DAD  D 

XCHG 


XCHG 

LHLD 

F 

DAD 

D 

MOV 

A,M 

I  NX 

H 

MOV 

H,M 

MOV 

L,A 

XCHG 

LHLD 

IN 

CALL 

$M9 

XCHG 

LHLD 

TO  2 

MOV 

A,M 

I  NX 

H 

MOV 

H,M 

MOV 

L,A 

DAD 

D 

DIVIDE  BY  128 

RL 

L 

RL 

H 

LD 

L,H 

LD 

A,0 

SBC 

A,0 

LD 

H,A 

SAVE  VALUE  OF 

X(I) 

.8080 

LDA  I 

;C  t’OR  MODEl:  X(1),CMP  0;  FOR  MODE  2:  X{2),CMP 


CMP 

0 

JZ 

DWN 

.Z80 

LD 

(XSAVE) ,HL 

.8080 

DWN: 

XCHG 

LHLD 

TOl 

MOV 

M,E 

I  NX 

H 

MOV 

M,D 

;F 

250  OUT=OUT- 

-G(I+1)*X(I+1) 

LHLD 

I 

DAD 

H 

SHLD 

XCHG 

TOO 

LHLD 

G 

DAD 

D 

SHLD 

TOl 

LHLD 

XCHG 

TOO 

LHLD 

X 

DAD 

D 

MOV 

A,M 

I  NX 

H 

MOV 

H,M 

MOV 

L,A 

72 


o  o 


8080 

V 


300 


Y(N)  =  XSAVE  = 


LOOP  N 


INC 

A 

PUSH 

HL 

PUSH 

BC 

POP 

HL 

LD 

c, (hl: 

ADD 

A,C 

LD 

C,A 

POP 

HL 

LD 

A,  (hl: 

OUT 

(C)  ,A 

DEC 

C 

INC 

HL 

LD 

A,  (hl: 

OUT 

(C)  ,A 

=  IN 

LHLD 

N 

DAD 

H 

DCX 

H 

DCX 

H 

XCHG 

LHLD 

Y 

DAD 

D 

SHLD 

TOO 

LHLD 

IN  ( 

=  X(l) 

OR  X(2) 

LHLD 

XSAVE 

XCHG 

LHLD 

TOO 

MOV 

M,E 

I  NX 

H 

MOV 

M,D 

LHLD 

N 

I  NX 

H 

I  NX 

H 

Y(N)=IN) 


C  LOAD  -(THE 


NUMBER 

LXI 

MOV 

MOV 

DAD 

XCHG 

JNC 


OF  TIMES  TO  LOOP(IOOO)  IN  HEX) 
B,-03E8 

D, H 

E, L 
B 

DON 


C  LOOP  K  REMOVE  SEMI-COLONS  FOR  LONGER  RUN  TIME 


C  LOAD  -(THE 


LHLD 
I  NX 

NUMBER  OF 


K 

H 

TIMES 


TO  LOOP(254)  IN  HEX) 


RETURN 


LXI 

B,-0 

MOV 

D,H 

MOV 

E,L 

DAD 

B 

XCHG 

JNC 

DOK 

74 


DSEG 


F: 

DW 

0 

G: 

DW 

0 

IH: 

DW 

0 

HX; 

DW 

0 

X: 

DW 

0 

Y: 

DW 

0 

I: 

DW 

0 

J: 

DW 

0 

K: 

DW 

0 

N: 

DW 

0 

OUT: 

DW 

0 

TOO: 

DW 

0 

VIN: 

DW 

0 

IN: 

DW 

0 

TOl: 

DW 

0 

T02: 

DW 

0 

VOUT: 

DW 

0 

FIFTY: 

DW 

50 

FORTY 8 : 

DW 

48 

ZERO: 

DW 

0 

VALUE: 

DW 

0 

CHAN: 

DW 

0 

BASE: 

DW 

0 

XSAVE: 

DW 

0 

END 

SUBROUTINE  ADIO 


.Z80 


NRDY: 


VALUE; 

CHAN: 

BASE: 

.ZBO 

ENTRY  DA 


DA: 


ENTRY  AD 

A/D  SERVICE  ROUTINE  -  CAPT  BRIGGS ( AFIT/ENY ,  WPAFB,OH) 
FORTRAN  CALLABLE:  CALL  AD ( VALUE , CHAN , BASE) 

GET  ONE  SAMPLE  FROM  THE  CHAN ' TH  CHANNEL 


ON  THE 

A/D  BOARD  WITH  BASE  ADDRESS  'BASE' 

LD 

(VALUE) 

,HL 

LD 

(CHAN)  , 

DE 

LD 

(BASE)  , 

bC 

EX 

DE  ,  HL 

; HL->CHAN 

LD 

A, (HL) 

,GET  CHAN  NO. 

LD 

HL,(BASE) 

LD 

C,(HL) 

;GET  BASE  I/O  ADDRESS  TO  C  REG  FOR  OUTING 

OUT 

(C)  ,A 

?MODE  0  TO  CHAN  NO. 

INC 

C 

;USES  BASE  ADDRESS  IN  C  REG 
; POINT  TO  START  CONVERSION  PORT 

LD 

A,0 

OUT 

(C)  ,A 

; START  CONVERSION 

DEC 

C 

; POINT  TO  BASE  REGISTER 

IN 

A,  (C) 

;GET  STATUS 

AND 

080H 

;BIT  7  IS  STATUS,  =1  IS  BUSY 

JR 

NZI,NRDY 

;NOT  ALL  O'S  =>  BUSY 

INC 

C 

; POINT  TO  BASE  ADD+1 

INC 

C 

; POINT  TO  DRL 

IN 

A,(C) 

;LOW  BYTE  OF  VALUE 

LD 

E,A 

INC 

C 

; POINT  TO  DRH 

IN 

A,(C) 

JHIGH  BYTE  OF  VALUE 

AND 

OFH 

;MASK  OUT  HIGH  NIBBLE 

LD 

D,A 

?DE=VALUE 

LD 

HL, (VALUE)  ;HL->WHERE  TO  PUT  VALUE 

LD 

(HL) ,E 

;PUT  LOW  BYTE  OF  VALUE 

INC 

HL 

LD 

(HL)  ,D 

;THAT  GIVES  THE  CALLER  THE  VALUE 

RET 

DVv 

0 

; STORAGE  FOR  ADDRESS  OF  VALUE 

DW 

0 

; STORAGE  FOR  ADDRESS  OF  CHANNEL  NO 

DW 

U 

; STORAGE  FOR  ADDRESS  OF  BASE  ADDRESS 

CALL 

DA(VAL, CHAN, BASE) 

CHAN 

IS 

0 

BASE 

IS 

72 (BASE  10) 

LD 

A,(DE) 

;GET  CHAN 

ADD 

A, A 

; DOUBLE  IT 

INC 

A 

;ADD  'jAE 

PUSH 

HL 

;SAVE  VAL 

PUSH 

BC 

• 

9 

POP 

HL 

;HL=>BASE 

LD 

C,(HL) 

;C=LOW  BYTE  OF  BASE 

LD 

C,A 

POP 

HL 

LD 

A,(HL) 

OUT 

(C)  ,A 

DEC 

C 

INC 

HL 

LD 

A,(HL) 

OUT 

RET 

(C)  ,A 

C=LOW  BYTE  VALUE  OF  PORT 

GET  VAL 

GET  LOW  BYTE 

PUT  LOW  BYTE 

C=HIGH  BYTE  PORT 

HL=>HI  BYTE 

GET  HI  BYTE 

PUT  HI  BYTE 


SUBROUTINE  PLOT 


SUBROUTINE  PLOT(Y) 

IMPLICIT  INTEGER(A-Z) 

DIMENSION  Y(IOOO) 

EXTERNAL  DA, WAIT 

C  PLOT  THE  VALUES  IN  ARRAY(Y(N)) 
YMX=0 

DO  200  N=l,399,2 
Y(N+1)=(Y(N+1) )/2 
IF(Y(N) .GT.O)  GO  TO  110 
YCOMP=-Y(N) 

GO  TO  115 
110  YCOMP=Y(N) 

115  CONTINUE 

IF  ( YCOMP.lt. YMX)  GO  TO  120 
YMX=YCOMP 

120  YMAX=2048+YMX+YMX/10 
YMIN=2048-YMX-YMX/10 
Y(N)=2048+Y(N) 
Y(N+1)=2048-Y(N+1) 

M=3 

C  PEN  UP 

CALL  DA(3500,3,72) 

C  WAIT  FOR  PEN  TO  COME  UP 
CALL  WAIT(M) 

C  GO  TO  Y-VALUE 

CALL  DA(Y(N) ,2,72) 

C  GO  TO  X- VALUE 

CALL  DA( Y(N+1) ,1,72) 

C  WAIT  FOR  PEN  TO  STOP  MOVING 
CALL  WAIT(M) 

C  PEN  DOWN 

CALL  DA(2048,3,72) 

CALL  WAIT(M) 

200  CONTINUE 

C  WHEN  DONE  PLOTTING , DRAW  BOUNDARIES 
C  GO  TO  THE  ORIGIN 

CALL  DA(3500,3,72) 

CALL  WAIT(M) 

CALL  DA( 2048 ,1 ,72) 

CALL  DA(2048,2,72) 

CALL  DA( 2048,3,72) 

CALL  WAIT(M) 

C  GO  AROUND  TWICE 
DO  220  1=1,2 
CALL  DA( YMAX,2,72) 

CALL  WAIT(M) 

CALL  DA(Y(N+1) ,1,72) 

CALL  WAIT(M) 

CALL  DA(YMIN,2,72) 

CALL  WAIT(M) 

CALL  DA(2048,1 ,72) 

220  CALL  WAIT(M) 


CALL  DA(2048,2r72) 
XM=2048 

V=(YMAX-2048)/30 

PV=2048+V 

MV=2048-V 

U=( YMAX-2048)/5 

DO  240  1=1,7 

C  INCREMENT  X-AXIS 
CALL  WAIT(M) 

CALL  DA( 2048 ,3 ,72) 

XM=XM+200 

CALL  WAIT(M) 

CALL  DA(XM,1,72) 

CALL  WAIT(M) 

CALL  DA(PV,2,72) 

CALL  WAIT(M) 

CALL  DA(MV,2,72) 

CALL  WAIT(M) 

240  CALL  DA(2048,2,72) 

C  RETURN  TO  ORIGIN 
CALL  WAIT(M) 

CALL  DA(Y(N+1) ,1,72) 
CALL  WAIT(M) 

CALL  DA(2048,1,72) 
CALL  WAIT(M) 

CALL  DA(2048,2,72) 
YM=2048-6*U 

C  INCREMENT  THE  Y-AXIS 
DO  270  1=1,10 
yM=YM+U 
CALL  WAIT(M) 

CALL  DA(YM,2,72) 

CALL  WAIT(M) 

CALL  DA(2055,1,72) 
CALL  WAIT(M) 

270  CALL  DA(2048 ,1 ,72) 
CALL  WAIT(M) 

CALL  DA(YMAX,2,72) 
CALL  WAIT(M) 

CALL  DA(2048,2,72) 
CALL  WAIT(M) 

CALL  DA(3500,3,72) 

RETURN 

END 


SUBROUTINE  VvRTBUF 


SUBROUTINE  WRTBUF ( BUF ,NS2 ,RUN ) 

IMPLICIT  INTEGER(A-Z) 

DIMENSION  BUF(l) 

DIMENSION  ZBUF(8) 

ST=0 

WRITE(1,50) 

FORMAT (3X, 'MODE  1 ' ,5X, 'TIME ' , 3 ( 5X , ' MODE  1 ' ,4X, 'TIME' ) ) 

N=NS2/8 

R=RUN+5 

DO  800  JA=1',N 

DO  90  C=l,8 

P=ST+C 

ZBUF(C)=BUF(P) 

WRITE (1,100)  (ZBUF(K)  ,K=1',8) 

WRITE(R,100) (ZBUF(K) ,K=1,8) 

FORMAT(/ ,4( 2X,I7 ,3X,I6) ) 

ST=ST+8 

CONTINUE 

NT=NS2-N*8 

IF(NT.EO.O)  GO  TO  910 
DO  890  D=1,NT 
P=ST+D 

ZBUF(D)=BUF(P) 

WRITE (1, 900 ) (ZBUF(M) ,M=1, NT) 

WRITE(R,900)  (ZBUF(M)i,M=l,NT) 

FORMAT(/,4(2X,I7,3X,I6) ) 

CONTINUE 

RETURN 

END 


SUBROUTINE  STCLK 

ENTRY  STCLK 

.Z80 

STCLK: 

LD  A,017H  CHANNEL  1  CTRL  WD  =TIME/16 
OUT  (079H),A 

LD  A, 09 AH  ;TIME  CONSTANT 
OUT  (079H),A 

LD  A,057H  , 'CHANNEL  2  CNTRL  WD=CTR 

OUT  (07AH) ,A 

LD  A,0FFH  ;TIME  CONSTANT* 2 56 ( BASE  10) 
OUT  (07AH),A 

LD  A,057H  ;CHANNEL  3  CNTRL  WD=CTR 
OUT  (07BH),A 

LD  A,0FFH  ;TIME  CONSTANT 

OUT  (07BH),A 

RET 

END 


I 


Turn  on  all  equipment  -  leave  High  Voltage  Dc  (HVDC)  off. 
Check  accelerometer  firmly  in  place. 

Check  all  connections  (especially  micro  dot)  secure. 
Assure  HVDC  dials  are  fully  CCW  and  turn  on. 

Select  operate  on  analog  computer/check  potentiometers. 
Simultaneously  set  both  HVDC  dials  at  200  ma. 

Turn  plotter  on,  center/check. 

Set  beam  initial  displacement  or  white  noise  generator, 
lace  disk  #3(operating  system)  in  drive  A  (label  faces 
right) . 

Place  disk  #5  (program)  in  drive  B  (label  faces  right). 
"Boot  up"  the  system  by  pressing  reset  button. 

Type;  B;  (return). 

Type;  TASK91  (return)  tor  mode  1  output,  and  TASK92  for 
mode  2  output.  TASK91  is  just  TASK9  linked  with  CNTL2M 
which  outputs  the  estimate  of  mode  1.  TASK92  is  formed 
similarly  using  CNTL2M  which  outputs  the  estimate  of  mode 
2.  Both  give  discrete  time  response.  For  the  continuous 
case  use  TASK91N  or  TASK92N  which  are  formed  by  deleting 
the  appropriate  ;  in  CNTL2M. 

Enter  run  number,  push  return,  and  release  the  initial 
displacement  mechanism,  or  push  RUN  tor  noise  input. 

Data  will  appear  on  screen  and  in  a  disk  file,  Fort0(run 
number  plus  tive).Dat,  as  tour  pairs  of  value  and  time 
(1  clock  cycle  =  .00005013  sec.). 


Appendix  7 

LiJbracy 

This  file  contains  a  list  of  the  information  stored  in 
room  147  pertaining  to  Dr.  Calico's  Flexible  Beam  Experiment. 
Vol.  I. 

A.  Microsoft  Fortran-80  Documentation  -  Microsoft. 

1.  Fortran-80  Reference  Manual  -  statements, 
functions,  syntax. 

2.  Fortran-80  User's  Manual  -  compiler  commands 
and  error  messages. 

3.  Microsoft  Utility  Software  Manual  -  Macro-80 
assembler*,  LINK-80  Linker  Loader , LIB-80  Lib 
Manager. 

B.  MC  5.16  Vibrations  Laboratory  Experiments  - 
Winter  1982. 

1.  Oscilloscope,  RMS  meter.  Oscillator. 

2.  Shaker,  Power  Supply,  Power  Amplifier,  Freq 
Counter. 

3.  Accelerometers,  Accelerometer  Amplifiers. 

4.  Unknown  Accelerometer  Calibration. 

5.  Cantilever  Beam  Response. 

6.  Measure  Resistance,  Output  Impedance. 

7.  Oscillograph,  Galvanometer  Operation. 

8.  Vibration  Decay  Experiment. 

C.  CP/M  Seminar  Lecture  Notes. 

Vol  II  EE  6.87  Microcomputer  Lab  -  Alan  Ross 

A.  Course  Notes  and  Assignments. 
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B.  ZCPR  Documentation. 

C.  SID  Symbolic  Intruction  Debugger  User's  Guide. 

D.  ZSID  Command  Summary. 

E.  MACRO-80  Assembler  Reference  Manual. 

F.  WM  Screen-Oriented  Text  Editor  User's  Manual 

E.  Manual  0;  S-100  Bus. 

F.  Manual  1:  Z80  CPU  Technical  Manual. 

G.  Manual  2;  Z80  CPU  Technical  Manual  and  Applica¬ 
tions  Note  (8251  Specification  Sheets  contained 
within) . 

H.  Manual  3:  Z80  PIO  Technical  Manual. 

I.  Manual  4;  Z80  CTC  Product  Specification. 

J.  Manual  5;  Z80  Family  Interrupt  Structure 
Tutorial. 

K.  Manual  6:  SBC-100  Single  Board  Computer. 

Vol  III.  Mainframe  Hardware. 

A.  Model  910  Terminal  -  Televideo. 

B.  S-100  Bus  Mainframe  -  Electronic  Control 
Technology . 

1.  Expandoram  II  -  Random  Access  Memory 
SD  Systems. 

2.  SD  Monitor  -Intruct ional  Publication. 

3.  SBC-100  Single  Board  Computer. 

4.  Mostek  3882  CTC  Technical  Manual  -  Mostek, 

5.  S-100  Bus  Extender  Board/Logic  Probe  Kit. 

6.  S-100  D/A  Converter  Board  -  Techmar  Inc. 

7.  S-100  A/D  Converter  Board  -  Techmar  Inc. 

8.  AIM-  12  Analog  Input  Module  -  Dual  Systems 
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Control . 

9.  AOM-  12  S-100  D/A  Converter  -  Dual  Systems 
Control. 

10.  HP3721A  Correlator  Documentation. 

11.  HP3720A  Spectrum  Display  Documentation. 

Vol  IV.  Background  Information. 

A.  PACOSS  -  Statement  of  Work  -  Passive  and  Active 
Control  of  Space  Structures. 

B.  DTIC  (Defense  Technical  Information  Center) 
Literature  Search. 

C.  ASIAC  (Aerospace  Structures  Information  and 
Analysis  Center)  Literature  Search. 

D.  User's  Guide  to  DTIC. 

E.  Control  of  a  Flexible  Satellite  Via  Elimination 
of  Observation  Spillover  -  Janiszewski/Calico  - 
AFIT. 

F.  Hardware  Demonstration  of  Flexible  Beam  Control 
Schaechter  JPL,  California  Institute  of 
Technology. 

G.  An  Experimental  Investigation  of  Modern  Modal 
Control-Herrick,  Canavin,  Strunce  -  Draper  Labs. 

Vol  V.  Tarbell  Disk  Control/Digital  Research  CP/M  Pubs. 

A.  Tarbell  Basic  -  Tarbell  Electronics. 

B.  Tarbell  Basic  15.4  Improvements  and  Changes, 

C.  Tarbell  CP/M  2.0  ( BIOS )/User ' s  Guide. 

D.  Tarbell  Double  Density  Floppy  Disk  Interface- 
Owner's  Manual. 

E.  Intro  to  CP/M  Features  and  Facilities-  Digital 
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F.  ED;  Context  Editor  for  theCP/M  DiskSystem 
User's  Manual, 

G.  CP/M  Dynaiti’ ::  Debugging  Tool  (DDT)-  User's  Guide. 

H.  CP/M  Assembler  (ASM)  -  User's  Guide. 

I.  CP/M  2  User's  Guide. 

J.  CP/M  2  Alteration  Guide. 

K.  CP/M  2  Interface  Guide. 

Vol  VI  Additional  Sof tware/Studor ' s  Thesis 

A.  TOTAL  USERS  MANUAL  -  An  interactive  computer- 
aided  design  program  tor  digital  and  continuous 
control  system  analysis  and  synthesis. 

B.  OPTCON  -  Solution  of  Linear  Quadratic  Optimal 
Control  Problem  -  Witt,  Mirmak,  &  Ray  March  1982 

C.  Thesis 

D.  Disk  Programs 

1.  Used  in  thesis 


2 


Future  applications 


George  Frederick  Studor,  Jr.  was  born  in  Syracuse,  New 
York  on  November  22,  1954.  Soon  after  high  school  graduation 
in  1972  he  began  his  Air  Force  service  as  a  cadet  at  the  USAF 
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in  August,  1977.  He  was  assigned  to  the  41st  TAS,  Pope  Air 
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